Lagrangian dynamics of thermal tracer particles 
in Navier-Stokes fluids 
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A key issue in fluid dynamics is the definition of the phase-space Lagrangian dynamics character- 
izing prescribed ideal fluids (i.e., continua), which is related to the dynamics of so-called ideal tracer 
particles moving in the same fluids. These are by definition particles of infinitesimal size which do 
not produce significant perturbations of the fluid fields and do not interact among themselves. In 
this work we point out that the phase-space Lagrangian description of incompressible Navier-Stokes 
fluids can be achieved by means of a particular subset of ideal tracer particles, denoted as thermal 
particles. For these particles the magnitude of their relative velocities - with respect to the local 
fluid velocity - is solely determined by the kinetic pressure, in turn, uniquely related to the fluid 
pressure. The dynamics of thermal tracer particles is shown to generate the time-evolution of the 
fluid fields by means of a suitable statistical model. The result is reached introducing a 1-D sta- 
tistical description of the fluid exclusively based on the ensemble of thermal tracer particles. In 
particular, it is proven that the statistic of thermal particles can be uniquely defined requiring that 
for these particles the directions of their initial relative velocities are defined by a suitable family of 
random coplanar unit vectors. 

PACS numbers: 51.50+v, 52.20-j, 52.27.Gr 



1 - INTRODUCTION 

The phase-space Lagrangian description of incompress- 
ible Navier-Stokes (NS) fluids is an essential ingredient 
of turbulence theory [l], 0]. Turbulence, in fact, can 
be viewed as the manifestation of stochastic behavior of 
the underlying phase-space dynamical system which ad- 
vances in time the state of the fluid (the so-called NS dy- 
namical system [3|). It is therefore important to provide 
an insight into the phase-space dynamics which charac- 
terizes these systems 0, IE El • 

In the framework of the theory of continua, for a given 
fluid system, there is a unique minimal class of such dy- 
namical systems which provide - via the introduction of 
an appropriate phase-space probability density function 
(PDF) - the time evolution of the complete set of fluid 
fields describing the same fluid. In turbulence theory this 
phase space is usually identified with the set T = £1 x U, 
with dim(T) = 6, denoted as restricted phase-space (with 
U = R 3 denoting the velocity space and r = O x U the 
closure of T ). Therefore it is natural to seek possible 
phase-space representations of this type for fluid systems. 

Although phase-space descriptions of incompressible 
NS fluids have been around for a long time, starting from 
the historical work of Hopf (Hopf, 1952 0]), Edwards 
(Edwards, 1964 0) and Rosen (Rosen, 1971 @]), a fun- 
damental issue concerns the construction of phase-space 
approaches exclusively based on classical statistical me- 
chanics (CSM). This concerns, in particular, the search 
of so-called inverse kinetic theories (IKT) able to yield 



the complete set of fluid fields describing the fluid. In 
these approaches, denoted as complete IKT's [13], the 
fluid is described in terms of a statistical ensemble of 
point-particles, called ideal tracer particles, whose phase- 
space trajectories are defined by the general solution of 
the initial-value problem of the form 



|r(t)=v(t), 



di 



£v(t)=F(r(t),t;/), 
r(i G ) = r D , 
v(* ) = v , 



(1) 



where F is a suitable vector field {mean-field force) and 
the state- vector x = (r, v) spans by assumption the phase 
space r = Q, x U. For them, as pointed out elsewhere [aja] 
and unlike traditional approaches 11, [13, 0, 14, H3], 
the dynamics (denned in terms of F) can be in principle 
uniquely established based exclusively on first principles. 
As a result, by construction, in a complete IKT, ideal 
tracer particles: 

• uniquely determine, by means of the 1-point statis- 
tics, the time evolution of the fluid fields; 

• interact with the fluid via a suitable mean-field 
force F; 

• do not perturb the fluid fields; 

• are collisionless (i.e., they do not interact among 
themselves); 
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• are subject to a mean-field force defined in such a 
way that the dynamics of the ensemble of tracer 
particles reproduces exactly, via a suitable statis- 
tics, the dynamics of the fluid fields. In particu- 
lar, the complete set of fluid fields is determined in 
terms of the 1-point PDF. 

An IKT of this type, in which tracer particles can have 
arbitrary initial velocities (i.e., which belong to the whole 
velocity space U = M 3 ), has been recently established for 
incompressible thermofiuids (see in particular, Tessarotto 
et al., [q]). In particular, in the case of an incompressible 
and isothermal NS fluid it is immediate to show that 
the dynamics of tracer particles generally differs from 
that of the corresponding fluid elements BE Eft The 
difference is intrinsic, i.e., it is produced by the mean- field 
force acting on the tracer particles (which differs from the 
local force density acting on the fluid elements), so that 
generally the fluid velocity is not a solution of the NS 
dynamical system. 

A peculiar characteristics of a complete IKT is, how- 
ever, that the mean-field force F contains an arbitrary 
scalar parameter (a) [13] • This feature, which - never- 
theless - is not reflected by the fluid fields (which remain 
by construction independent of the same parameter), can 
also be interpreted as the non-uniqueness in the defini- 
tion of the NS dynamical system. 

An interesting question is, however, whether a proper 
subset of solutions (of the NS dynamical system) exists 
which can be used to characterize uniquely incompress- 
ible NS fluids and resolve also uniquely the previous inde- 
terminacy issue. This refers, in particular, to the possi- 
bility of identifying a proper subset of the family of tracer 
particles, characterized by a phase-space of lower dimen- 
sion (i.e., smaller than 6) and whose states uniquely de- 
termine exactly - via a suitable statistics - both the com- 
plete set of fluid fields and their time evolution. 



Goals of the investigation 

Purpose of this paper is to point out the discovery 
of a family of tracer particles denoted as thermal tracer 
particles, which have all such features. These tracer par- 
ticles are characterized by relative velocities belonging 
to a subset of the velocity space U of dimension I. As 
a main consequence, it is shown (see THM.f) that the 
dynamics of these particles uniquely determines - via a 
suitable statistics - the time-evolution of the NS fluid. 
The ensemble of these particles is defined by the require- 
ment that for each tracer particle the magnitude of its 
relative velocity (with respect to the local fluid velocity) 
is proportional to the fluid pressure. More precisely, de- 
noting by x(t) = {r(i),v(i)} the phase-space trajectory 
of a generic thermal tracer particle, with r = r(t) and 



v = v(t) = 3jr(t) its position and velocity, V (r,t) the 
local fluid velocity and u = v V(r, t) the relative parti- 
cle velocity, we intend to prove that *-th(t) = {r(t), v(t)}, 
with 



where 



5 r(*) = v(t), 
v(f)=V(r,t)+u(t), 



u(i) = n(r,i)u(t), 
u(t) = v th (v,t) = 



'2pi(r,t) 



Po 



(2) 
(3) 

(4) 
(5) 



is a particular solution of the NS dynamical system. In 
particular, v t h denotes the thermal velocity produced by 
the kinetic pressure pi(r, t) [which is, in turn, related 
to the fluid pressure, see Eq. (f24[) below], while n(r,t) is 
the rotating unit vector which satisfies the initial-value 
problem 



£n(ri,t) = n(r,t) x fi(r,i), 
n(r(t G ),t D ) = n(r ,t ). 



Here fi(r, t) is the angular velocity 
fi(r,t)=a£, 



(6) 



(7) 



produced only by the fluid vorticity £ =V x V, while the 
unit vector n(r, t) is subject to the requirement 



n(r,*)-Vp 1 (r,t) = 



(8) 



{constraint of orthogonality) . In Eq.([7| a = a(r, t) is the 
real scal ar p arameter which characterizes the mean-field 
force F [17]. The arbitrariness of a implies that the def- 
inition of the angular velocity S7(r, t) is non-unique, a 
feature which can be viewed as as kind of kinematic in- 
determinacy of the theory. In Ref. [17j a was determined 
based on phenomenological arguments as 



a = 1/2. 



(9) 



Here we intend to show, however, that the requirement 
of the existence of the solutions of the type ©-J5]) ac- 
tually imposes also a kinematic constraint on the same 
parameter, i.e., that there results also 



r, d 
n - V dt Pl 



0, 



(10) 



where uj = (V x V) x Vpi- In other words (see THM.l) 
the parameter a can also be defined in such a way that, if 
at the initial time t , n(r ,t ) satisfies the constraint of 
orthogonality (with respect to V j>i(r , i ), then at any 
time t, n(r(i),i) will satisfy the same constraint provided 
(JTUJ) is satisfied too. Therefore, in such a case the rela- 
tive velocity of a generic thermal tracer particle u(i) is 
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a rotating vector orthogonal to Vf>i(r, t). In particular, 
it is found that while its magnitude u(t) depends, ac- 
cording to Eq.©, only on the kinetic pressure pi(r,t), 
the angular rotation velocity i~2(r, t), which characterizes 
the direction of the relative velocity, is produced by fluid 
vorticity. 

Main consequence of the theory here presented is that 
dynamics of an incompressible NS fluid, i.e., the fluid 
fields and the complete set of fluid equations, can be rep- 
resented solely in terms of the kinetic state of the ther- 
mal tracer particles (see THM.2). The goal is reached 
by proving that the dynamics of these particles uniquely 
generates, via a suitable 1-D statistics, the time-evolution 
of the fluid fields. 



Scheme of presentation 

The general framework (of this work) is provided by 
the phase-space representation for an incompressible NS 
fluid previously developed (see in particular Tessarotto 
and coworkers [R IR HR flR] ) . In Sec. 2 the relevant aspect 
of the theory are recalled, with particular reference to the 
so-called Hopf- Rosen- Edwards (HRE) approach @, H, 
(see also Monin and Yaglom, 1975 [l9| and Pope, 2000 
(H) and the IKT approach developed in Refs.(3,[l3|. This 
permits us to identify the NS dynamical system, via the 
construction of a suitable mean-field force. After dis- 
playing the corresponding phase-space Lagrangian for- 
mulation (in Sec. 3), the existence of the thermal tracer- 
particles solutions is proven in THM.l (see Section 4). 
Finally the statistics of thermal tracer particles is inves- 
tigated in Sec. 5 (THM.2). This is shown to depend 
only on the direction of the tracer-particle relative veloc- 
ity, defined by a coplanar family of unit vectors. 



2 - PHASE-SPACE APPROACHES TO FLUID 
DYNAMICS 

A fundamental aspect of fluid dynamics is the con- 
struction of phase-space approaches based on classical 
statistical mechanics (CSM). Indeed, phase-space tech- 
niques are well known both in classical and quantum fluid 
dynamics. In this connection, a particular viewpoint is 
represented by the class of so-called inverse problems, 
involving the search of an inverse kinetic theory (IKT) 
able to yield the complete set of fluid equations for the 
fluid fields, via the introduction of a suitable phase-space 
PDF /(x,t; Z) on the restricted phase-space V [l8]. This 
requires, by assumption, that the PDF must depend at 
least functionally, i.e., via appropriate "moments" , on the 
(deterministic) fluid fields {Z} characterizing the fluid 
(or more generally on a suitable subset of {Z}). 

A particular realization for /(x,i;.Z), holding for an 
incompressible NS fluid, is provided by the well-known 



Hopf- Rosen- Edwards (HRE) approach 0, d, (see also 
Monin and Yaglom, 1975 and Pope, 2000 Q). In this 
case /(x,i; Z) satisfies the Liouville equation [or inverse 
kinetic equation] 

L(r,v,t;ZJ)f(x,t;Z)=0, (11) 

with /(x,t; Z) = f H (x,t; Z) denoting the PDF 

f H (x,t;Z) = 5{v-V{v,t)) (12) 

and L(r,v,t; Z, f) a suitable Liouville streaming opera- 
tor. Here the notation is standard [3j. Thus x = (r, v) is 
a state- vector which spans the phase space T = 57 x U, 57 
and U = R 3 denoting respectively the configuration and 
velocity spaces. Moreover, L(r,v,t; Z, f) reads 

L(r, v, t; Z, /). = _.+_ . {X(x, t; Z)-} = (13) 

4-+v-V. + A. { F(x,,Z). } 

and X(x, t; Z) is the vector field 

X(x,<;Z) = {v,F}. 



(14) 



In particular, in the HRE approach F is identified with 
F = Fjj, i.e., the total fluid force per unit mass acting on 
each fluid element [see Eq. (f79")) in Appendix A] . It follows 
by construction that: 

• the first two velocity moments of /(x,i; Z) are nec- 
essarily 



J dvf(x,t;Z) = 1, 

u 

J dw/(x,i;Z) = V(r,i); 



(15) 
(16) 



• the corresponding moment equations obtained 
above coincide manifestly with the complete set of 
fluid equations [see Eqs. ([73| and (|74p]: 

• since the PDF f(x,t;Z) of Eq. (fl"2|) . is independent 
of the fluid pressure, p(r,t) cannot, manifestly, be 
represented as a moment of the same PDF. 

However, a more general viewpoint is represented by 
the class of so-called complete IKT's able to yield as 
moments of the PDF the whole set of fluid fields {Z} 
which appear in the fluid equations. A theory of this type 
has been recently developed by Tessarotto and coworkers 
[H,[l8|. In the case of INSE the theory must include neces- 
sarily, among the PDF moments, also the fluid pressure. 
This requires (Principle of correspondence) that, besides 
Eqs.(fT5|) and (|16|) . there must exist a suitable phase-space 
function G(x,t) for which it results identically in 57 x / 



p(r,t) = / dv/(x,t;Z)G(x,t). 



(17) 
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In particular, one can show that a possible choice is pro- 
vided by the position Q 



G(x,i) 



PoU 



p (t) - (/)(r,t), 



(18) 



where u is the relative velocity u = v — V(r,t) and 
Po(t) > is a strictly positive smooth real function to 
be defined (see below). 

The actual construction of a complete IKT of this 
type relies on classical statistical mechanics [13, l2fij |. 
and its formulation comprises the usual axioms of CSM. 
It follows that, again, /(x,i; Z) must fulfill necessar- 
ily the Liouville equation (which can be intended 
as a phase-space Eulerian equation). Here by con- 
struction X(x, i; Z) is so defined that the fluid fields 
{V (r ,t) , p(r ,t)} must be uniquely determined by the 
moment equations (|16|) and (I17j) . It is important to 
stress that, by suitably specifying the functional class 
{/(x,t; Z)}, the actual form of the vector field X(x, t; Z) 
can be explicitly constructed Q and proven to be unique 
US El . Thus, a particular solution of Eq. (fTTI) can be 
obtained by requiring that /(x,i; Z) for any (x,t) E V x I 
satisfies the two entropic principles: 

• Axiom of maximum entropy (also known as princi- 
ple of entropy maximization (PEM); Jaynes, 1957 
[22|), This requires that: 1) /(x,i; Z) be a strictly- 
positive ordinary function which admits for all t 6 I 
the moment 



S(f(t)) 



J dx/(x,t;Z) In/(x,t;Z), 



(19) 



S(f(t)) denoting the Boltzmann-Shannon (BS) en- 
tropy associated to the PDF f(x,t;Z); 2) that for 
all t e / the probability density /(i) = f(x,t,Z) 
satisfies the constrained maximal variational prin- 
ciple 



(20) 



SS(f(t)) = 0, 



with f(t) required to satisfy the constraints (jT5|) 
and (fT6|) and (flTf : 

• Axiom of conservation of entropy (also known as 
constant H-theorem 17|). This requires that there 
results identically for all t E I 



di 



S(f(t)) = 0. 



(21) 



Imposing PEM requires necessarily that /(x,i; Z) must 
coincide with a local Maxwellian distribution [3, Irjj 



/ M (u;pi(r,i)) = 



1 



ir 2 vth(r,t) 



exp 



ut/t(r,t) 



(22) 



where u 



V(r,t), 



wtfc(r,t) = y/2 Pl (r,t)/ Po , (23) 

is the thermal velocity due to the kinetic pressure 

Pl (r, t) = Po (t) + p{r, t) - 4>(r, t) , (24) 

and the pseudo-pressure Po(t ) > is uniquely deter- 
mined in such a way to assure that in the time interval / 
the constant H-theorem [Eq. (|2"Tj) ] is satisfied identically. 
As a further consequence, one can prove [l7| that the 
requirement that 



f( X ,t;Z) = f M (u;p 1 (v,t)) 



(25) 



be a solution of Eq. (fTTj) implies necessarily that the vec- 
tor field F depends functionally on fu and {Z} , hence 
it can be intended as a mean-field force [l^|. As pointed 
out previously [13], the form of F is actually non- unique. 
More precisely it can be cast in the form: 



F(x,£; /m; Z; a) = F (x,t; f M ; Z; a) 



¥ 1 {x,t;f M ;Z), 
(26) 



1 



F (x,f; / M ; a) = —f R + D(r, t; a)+//V 2 V, 
Po 



Fi(x,t;/ 



1 



Mj 



uA(r,t;/) + AFi(x,t;/. 



M) 



AFi(x,t; /m) = ^— Vpi <^ 

2pi L 2 

where D(r, t; a) is the convective term 

D(r, t;a) = u- VVfe + aVV • u, 



(27) 
(28) 
(29) 

(30) 



b = 1 — a, and a = a(r, t) is an arbitrary scalar real func- 
tion. Finally, the vector field A(r,t; /m) is now defined 
as 

A{T,t;f M ) = ~Pi- (31) 
Pi at 



Po 

Pl 



V 2 /2 + V • (VV 2 /2) V • f - z/V-V 2 V 



Po 



1 

Pl 



D 
Dt 



Pi 



PoV 2 



V- V 



where a term proportional to V • V [and hence identically 
vanishing when the isochoricity condition is satisfied; see 
Eq. ([73]) in Appendix A] has been added. 
Here we stress that: 



1. as previously pointed out 21( the convective term 
D(r, t;a) is non- unique due to the indeterminacy 
of the parameter a. This implies also the non- 
uniqueness of the NS dynamical system (see next 
Section). The choice (j9]) for the parameter a then 
corresponds to the symmetry condition [2lJ 



a = b = 1/2; 



(32) 
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2. the complete IKT holds in principle for arbitrary 
choices of a, to be considered either a real constant 
or a smooth real function a(r, t); 

3. as indicated elsewhere [3j the present theory can 
also be generalized in a straightforward way also 
to non-Maxwcllian initial distribution functions 
/(x,t;Z). 



3 - PHASE-SPACE LAGRANGIAN 
FORMULATION 

The results of the previous Section permit us to formu- 
late in a straightforward way the phase-space Lagrangian 
formulation for an incompressible NS fluid. The complete 
IKT implies, in fact, that there exists a phase-space clas- 
sical dynamical system (called NS dynamical system) 

x -» x(t) = {r(i), v(t)} = T Mo x ee x ( Xo , i OJ t), (33) 

Tt,t denoting the corresponding evolution operator, 
which generates the time-evolution of the complete set 
of fluid fields, i.e., is such that 



{Z((r(t),t)} = T ttto {Z (v )} 



(34) 



[here Z (r ) = Z(r ,t ) denotes the initial fluid fields 
evaluates at t = t Q ; see Eq. ([75|) ]. This means that the NS 
dynamical system is identified with the flow generated by 
the vector field X defined by Eq.([T3|l. In view of Eqs. (|26]) - 
(13T1) this must be considered of the form 



X(x,i;/;Z;a). 



(35) 



Hence, x(t) is obtained as the solution of the initial- 
value problem defined by Eq.([T]), where, in particular, 
the mean-field force F is expected to differ generally from 
Fh as defined by Eq.([79| [see Appendix]. Moreover the 
state- vector x = (r, v) spans by assumption the phase 
space r ee Q x U, while and U denote appropriate con- 
figuration and velocity spaces. The integral curves x(t) of 
Eq. {!]) can therefore be interpreted as Lagrangian phase- 
space trajectories of tracer particles, whose dynamics is 
determined uniquely by the time evolution of the fluid. 
The latter can also be cast in an equivalent phase-space 
Lagrangian form. By assumption the initial-value prob- 
lem ([T]) is well-posed, namely it defines a suitably smooth 
diffeomorphism. Here, by construction: 

• denoting by x(i) = x( x o,£o,£) the solution of the 
initial-value problem ([1]), x G = x(x(i),i, i ) is its 
inverse. Both are assumed to be suitably smooth 
functions of the relevant parameters; 

• both x(i) = x( x o, t , t) and x c = x(x(t), t, t Q ) iden- 
tify admissible LP's of the dynamical system; 

• r(t) is the Lagrangian trajectory which belongs to 
the fluid domain Q: 



• v(i) and F(r(t),t;f) are respectively the La- 
grangian velocity and acceleration, both spanning 
the vector space R 3 . In particular, F(r(t),t;/), 
which is defined by Eqs.([26 | - ((3T|) . and depends 
functionally on the kinetic probability density 
f(x,t), is the Lagrangian acceleration which cor- 
responds to an arbitrary kinetic probability density 
/(x,i); 

• /(x, t) is a particular solution of the inverse kinetic 
equation (fTT|) . 



3x(t) 
dx 



of 



It follows that the Jacobian [J(x(t),i) = 

the map x Q — ► x(t), which is generated by Eq.([T]) for 
/ ee /m(x,*) reads 



J(x(t),t) 

Vth(t) 
Vth{to) 



dx(t) 



exp 



dx 

u 2 {t ) | u 2 (t) 

Vth(to) V t h{t) 



(36) 



Furthermore, it is immediate to prove that the inverse 
kinetic equation can be written in Lagrangian form. This 
yields the so-called integral Liouville equation 



J(x{t),t)f(x(t),t;Z)=f{x ,t ,Z), 



(37) 



with /(x D , t a , Z) denoting the initial PDF and 

J(x(t),i) = the Jacobian of the map ([33]) . As 

a consequence, for all (x, t) s T x I (with T = ft x U) the 
PDF /(x, t; Z) can be represented explicitly in terms of 
the initial PDF and reads 



fMZ) 



i 



J(t;Z) 



f(x( x , t , t ), t o,Z), 



(38) 



where J(t;Z) ee J(x,t), X (x, t, t) ee T^ t x and x ee x(£). 
In particular, in terms of Eq. ((38|) the fluid fields can be 
explicitly evaluated, yielding: 



V(r,t) 



(ivv- 



-f( X (x,t,t),t ,Z), (39) 



J(t;Z)' 

P(r,t) = J dyj^—f( X (^t,t),t ,Z) (40) 



p (0-^(r,i) 



From the mathematical standpoint main consequences of 
the IKT representation are that: 1) the Lagrangian for- 
mulation (of IKT) is uniquely specified by the proper 
definition of a suitable family of phase-space LP's; 2) 
Eq. ([38]) uniquely specifies the time-evolution of the Eule- 
rian PDF, /(x(£), i), which is represented in terms of the 
initial distribution function / (x ) and the LP's defined 
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by the INSE dynamical system; 3) the time-evolution 
of the fluid fields {p = p > 0, V,p > 0,T > 0, S T } is 
uniquely specified via the PDF /(x(i), *); 4) Eq.® also 
provides the connection between Lagrangian and Eule- 
rian viewpoints. In fact the Eulerian PDF, /(x,t), is 
simply obtained from Eg. ([55]) by letting x = x(£) in the 
same equation. As a result, the Eulerian and Lagrangian 
formulations of IKT, and hence of the underlying moment 
(i.e., fluid) equations, are manifestly equivalent. From the 
physical viewpoint, it is worth mentioning that the LP's 
here defined can be interpreted as phase-space trajecto- 
ries of the particles of the fluid, to be considered as a 
set of "classical molecules", i.e., point particles with pre- 
scribed mass, which interact with the fluid only via the 
action of a suitable mean-field force kind. The ensem- 
ble motion of these particles has been defined in such a 
way that it uniquely determines the time evolution both 
of the kinetic distribution functions and of the relevant 
fluid fields which characterize the NS fluid. 



4 - THE DYNAMICS OF THERMAL TRACER 
PARTICLES 

An interesting issue is the search of possible exact par- 
ticular solutions of the Lagrangian equations |T]). The 
general solution of Eq. (JTJ) implies 

v(t) = V(r(t), t) + n(r(i), £>(r(i), £), (41) 

where n(r(i), t) is the rotating unit vector which satisfies 
the initial- value problem Eq.®. For arbitrary tracer- 
particle solutions it is obvious that n(r(i),£) is an arbi- 
trary unit vector belonging to the unit sphere [5(E/)] of 
velocity space U = K 3 . Hence, by definition, S(U) is a 
set of dimension 2. In fact, Eq.(JTJ) yields for the relative 
velocity u = v V(r, t) the equation 



— u = F u (x,t;f M ;Z;a), 



(42) 



F u (x,i; fu\ Z\ a) = aVV • u — au • VV+ (43) 



u 

Vth 

2pi 

It follows 
d u 2 u 2 1 

dt T~Tp7 



D 
Dt 



Pi 



u 

Vth 



. PoV 2 

2pi 
2 1 
2 



V- V 



PoV 2 



V • V 



1 

Wo 



u-Vpi. (44) 



df 1 ' 2 

Furthermore, Eg . ([42"]) . in view of Eq. (l44]) . implies for n 
d 



dt 



n(r£,t) = n(r,i) x (l(r,t)- 



(45) 



1 f u 2 1 

{ o \ n ( r ' *) X t n ( r ' *) X V ^i] 

Po I vth 2 



There it follows: 

Theorem 1 - Existence of thermal tracer parti- 
cles 

If the fluid fields {V(r, t),p(r, t)} are classical solu- 
tions of INSE [see Eqs. |7£| )- |76] ) in Appendix A], then 
a particular solution x(i) = x.Th(t) of Eq.fcfy is provided 
by Eqs.^-^, subject to the orthogonality and kinematic 
constraints 0) and H0\) . 

PROOF 

First we notice, invoking Eq. ([44|) and assuming that 
V(r, t) satisfies in particular the isochoricity condition 
[Eq. ([74]) , in Appendix A] , that Eq. implies necessarily 



d u 2 



1 d 



dt 2 2 pi dt 



Pi- 



(46) 



It follows that the unit vector n (r, t) must satisfy iden- 
tically Eq.©. On the other hand, thanks to Eq.©, ac- 
cording to Eq. ([45|) . n (r, t) is a rotating unit vector with 
angular velocity fi(r, t) defined by Eq.©- Hence, the 
condition of orthogonality Eq. © requires 



d d 
—n • Vpi + n- — Vpi = 

dt dt 



(47) 



too. Eq.© then implies manifestly also Eq.fTU]), which 
can be viewed as a kinematic constraint for the parameter 
a, to be considered as dynamical function. If uj is non- 
vanishing and the scalar product n ■ u> is defined so that it 
is 7^ it is obvious that the previous equation can always 
be satisfied. On the contrary, if locally u> = 0, this can 
happen either if 1) £ = 0, 2) Vpi = or 3) £ || Vpi- 

In the first case, it follows Eq.© reduces to = 
which requires also n • Vpi = and moreover V is po- 
tential, i.e., there exists a real scalar G(r,t) such that 
V = VG(r,£). This implies that G(r,t) must satisfy the 
PDE 



l V G + iv(|VG| 2 ) +Vp! = 0, 
namely locally it must be 



at 



VG 



-V 



(|VGf 







(48) 



(49) 



too. In the second case the direction of n is arbitrary 
so that it can always be chosen in such a way to satisfy 
Eq.©. Finally, in the third case if n is orthogonal Vpi it 
remains always so by construction thanks to Eq.©. Let 
us now impose the constraints ©, pH]) and consider an 
arbitrary initial condition of the form 



with 



x (*o) = {r(0 = r , v(t Q ) = v } 



v Q = V(r Q , t ) + n(r Q , t )u(t ), 



(50) 



(51) 



u(t ) = v t h{r ,t ) = \ , (52) 



Po 



n(r ,t ) ■ V pi(r ,i ) = 



(53) 
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Here {V(r D , t ),p(r , t a )} are arbitrary initial fluid fields 
consistent with INSE (see Appendix A), while pi(r 0) i ) 
and V pi(r ,t ) are defined according to Eqs. (f2"4")) . Fi- 
nally, n(r ,t ) = n is an arbitrary unit vector orthog- 
onal to V pi(r ,i ). It is obvious that the particular 
solution x(i) which satisfies the previous initial condi- 
tions (|50|) - (f5"3|) coincides necessarily with xtt(^) defined 
by Eqs.Q-d). Q.E.D. 

Let us analyze the physical interpretation of the theo- 
rem. 

A basic consequence is that an arbitrary tracer parti- 
cle, subject to the kinematic constraint (|10|) and fulfilling 
the initial conditions (|50|) - ((53|) is necessarily a thermal 
tracer particle for all £ /. Hence, at time t its relative 
velocity is 



u(r(t ),t)=n(r(t ),t)* 



'2 Pl (r(t),t) 



Po 



(54) 



with n denoting a rotating unit vector in accordance to 
Eqs.f6|) and j7]). The consequence is that tracer particles 
with initial velocity (|50I) - (I53I) behave as thermal parti- 
cles, i.e., at any time t their relative velocity [with re- 
spect to the local fluid element] is solely determined by 
the kinetic pressure p\ (r, t) [in turn related via Eq. (|24]) 
to the fluid pressure]. The direction [n(r, t)] of their rel- 
ative velocity u(r, t) depends, instead, on the angular 
velocity Jl(r, t) defined by Eq.((7j) and depending only on 
the fluid vorticity. However, in order to assure the con- 
dition of orthogonality ([5]) for the relative velocity, the 
kinematic constraint (fTU|) must be satisfies so that both 
f2(r, t) and the mean- field force F are both uniquely de- 
termined. This implies that u necessarily belongs to the 
I-dimensional subset of velocity space U t h C U for which 
the orthogonality condition Eq.© holds. 



THERMAL TRACER PARTICLE 1-D 
STATISTICS 



Apart the physical insight afforded by the qualitative 
behavior of tracer-particle dynamics, the existence of the 
particular solution defined by Eqs.©-© is important 
in order to establish a phase-space statistical description 
for incompressible NS fluids. In this section we intend 
to prove that the dynamics of thermal tracer particles is 
actually sufficient to determine uniquely the time evolu- 
tion of the fluid fields. In other words, to characterize 
completely the statistical behavior of the fluid it is suf- 
ficient, in principle, to determine only the dynamics of 
these tracer particles. 

In view of the general form of the solution for ther- 
mal tracer particles [see Eqs.([2])-([7])], it is obvious that 
the statistics depends only on the unit vector n(r, t) 
which defines the direction of the particle relative ve- 
locity with respect to the fluid [see Eq.©]. Thanks to 



the orthogonality constraint ([8]) placed on the dynam- 
ics of thermal tracer particles it follows that n(r, t) is 
orthogonal to Vpi(r,i) = e z \\/p%(r, t)| , i.e., n(r,i) must 
belong to a family of coplanar unit vectors which are tan- 
gent to the isobaric surface {pi(r,£) — const.}. Hence, 
the only admissible velocity-space statistics is necessar- 
ily one-dimensional (since it is defined on the subset of 
velocity space U t h define above). This is obtained by 
considering n(r, t) as a stochastic unit vector character- 
ized by a suitable stochastic PDF g. It is obvious that 
g is necessarily related to the PDF introduced in the 
complete IKT described in the previous section, i.e., to 
/(x, t; Z). In case (|25j) . by parametrizing n in cylindrical 
coordinates 



n = n (r,$, t) = e x cos $ + e y sin 



(55) 



it follows that g must be necessarily identified with the 
uniform distribution 

g = 1/2tt. (56) 
Thus, introducing the stochastic average 



r 27T 

= / d»g 



(57) 



[where the angle-integration is performed at constant 
r = r(t ) and t] it follows by definition that 

(I) = I, (58) 
(n) = 0, (59) 

(nn) = -| 2 = - [e x e x + e y e y ] . (60) 
As a consequence it is immediate to show that: 

• the fluid fields {V(r, t),p(r, t)} can be related to 
suitable stochastic averages of thermal tracer par- 
ticle dynamics; 

• similarly, by taking the stochastic average of the 
Lagrangian equations (fTJ) and (|4"4"1) , the stochastic- 
averaged equations 



-<v(t)> = <F(x,i;./M,Z)}, 



(61) 



™(« 2 (*)) = <u(£).F u (x,i;/ M ,Z)> (62) 

can be evaluated explicitly for thermal tracer par- 
ticles. 



The following result holds. 

Theorem 2 - 1-D statistics of thermal tracer 
particles 

In validity of THM.l, let us assume that x(f) = 
{r(t),v(t)} is an arbitrary phase-space trajectory of a 
thermal tracer particle [i.e., a solution of the form given 
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by Eqs. QJ)-0)/. TTien, due to fH and (5?|) ; it follows 
that: 

1) the stochastic averages (v(*)) and (u(i) 2 ) are 
related to the local values of the NS fluid fields 
{V(r, t),p(r, t)} , namely it results 



<v(t)> = V(r,i), (63) 



^(u(t) 2 )-po(t) + <f>(r,t) = p(r,t), (64) 

where po(t) and <f>(r,t) are defined by J£^ ),(fff7j )) [see Ap- 
pendix A] and the constant H-theorem i21\) ; 

2) equations 161]) and i62\) deliver respectively the NS 
equation and the isochoricity condition [see Eqs. FW an d 
Eq. (73\) in Appendix A ]. 

PROOF 

The proof of Eas. l|63|) and (|64[) follows immediately 
from Eqs. j58|)- (J60j) and Eqs. 

To evaluate the stochastic average of Eq. (f6l"j) we first 
recall the identity 

d d d 
— v=— V+ (V + mt)-VV + — u. (65) 

Lib L/6 (JjO 

Invoking Eqs.gU) and |g3J) it follows 



iv) = |v(r,t) + V.VV + (66) 



u) = Fh + (F u (x,t; /m; •Z')) , 



where F# is given by Eq. ([79|) [see Appendix A] and 

(F u (x,t;f M ; Z)> = ^V P1 {— - \ ) = 0. (67) 
2pi [v t h 2 J 

Hence Eq. ([66|) [and Eq. lfBTj) ] coincides manifestly with 
the NS equation (f74]l . 

In a similar way it is immediate to evaluate Eq. (|62|) . 
From Eqs.(U]) and (O, in validity of the constraint (jHJ), 
one obtains in fact: 



d u 2 

j t - ) =(u-F u ( X ,t;f M ;Z)), (OM 



while under the same conditions from Eq. (I43|) it follows 
*!r)-(..!v) + «V + o>.W. + 

= — ^v-v+ — — 

2 p -Dip 



d u 2 \ _ D pi(r,i) 



Eqs. (|62|) and (|68|) necessarily yield the identity 

— v-v = o. 



(70) 



(u-F„(x,i;/ M ;Z)) = — 



V- V 



(69) 



Hence, provided V 2 ^ 0, also the isochoricity condi- 
tion [see Eq. (fT3"l) in Appendix A] is identically fulfilled. 
Q.E.D. 



6 - CONCLUDING REMARKS 

In this paper properties of the NS dynamical system, 
advancing in time the state of an incompressible NS fluid, 
have been investigated. This refers, in particular, to the 
proof of the existence of a particular subset of phase- 
space trajectories (solutions of the NS dynamical) repre- 
senting the dynamics of so-called thermal tracer particles. 
The result has been reached by imposing suitable con- 
straints to the direction of u (relative velocity) of tracer 
particles and on the parameter a which characterizes its 
time evolution. In turn, the same constraint determines 
uniquely also the mean-field force F which defines the 
NS dynamical system. As a main consequence, we have 
proven (THM.l) that there exists a family of particu- 
lar solutions x(t) = {r(t),v(t)} which have the following 
features: 

1. the magnitude of the relative velocity with respect 
to the fluid, u = mt, is determined only by the 
kinetic pressure p\ (i.e., tracer particles of this type 
behave as thermal particles); 

2. the direction of the relative velocity n is a rotating 
unit vector. Its angular velocity is due to the fluid 
vorticity [see Eq.0]; 

3. the dimension of the subset of the velocity space 
Uth Q U spanned by these solutions is dim(Uth) = 
1; hence their phase space T t h—^ x Uth is of dimen- 
sion dim(r t ^) = 4. 

Furthermore, the TD velocity-space statistics of ther- 
mal tracer particles has been investigated (Section 5). As 
a main consequence (see THM.2) : 

1. The state of thermal tracer particles has been 
proven to determine uniquely, via a suitable statis- 
tics [namely by means of the stochastic averages 
(|63p and ([64]) ]. the complete set of fluid fields 
{V,p} ■ As a consequence the state of the fluid has 
been proven to be uniquely related to that of the 
thermal tracer particles. 

2. The complete set of fluid equations has been repre- 
sented in terms of stochastic averages of the phase- 
space Lagrangian equations for thermal tracer par- 
ticles (ED and 
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These conclusions show that a phase-space Lagrangian 
dynamics for an incompressible NS fluid based exclusively 
on the dynamics of thermal tracer particles is possible. 
The result appears relevant because of its outstanding 
physical interpretation. The phase-space Lagrangian dy- 
namics here determined permits, in fact, to advance in 
time self-consistently the fluid fields, i.e., in such a way 
that they satisfy identically the required set of fluid equa- 
tions. The proof of this statement is straightforward and 
is given in Appendix B. 
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APPENDIX A - DETERMINISTIC 
DESCRIPTION OF NS FLUIDS 

In fluid dynamics the state of a fluid is assumed to 
be prescribed by an appropriate set of suitably smooth 
functions {Z} = {Zi,i= l,n} denoted as fluid fields. 
These are required to be real functions of the form 



V(r, t) and p(r, t) > denoting respectively the deter- 
ministic fluid velocity and pressure and by the incom- 
pressible NS equations (INSE): 



P 

V- V 
AV 

Z{r,t ) 
Z(r,t)\ dn 



Po, 

0, 

o, 

Zo(r), 

Zw(Y-,t)\ da . 



Here N is the NS nonlinear operator 



D 



' H 



1 

Po 



-Vp- 



1 

Po 



-f + uV 2 V, 



(72) 
(73) 
(74) 
(75) 
(76) 



(77) 
(78) 
(79) 



with j^;V,F H and f(r,t) denoting respectively the La- 
grangian fluid acceleration, the total force per unit mass 
and the volume force density acting on the fluid, while 
p > 0, v > the constant mass density and the con- 
stant kinematic viscosity. In particular, f is assumed of 
the form 



f = -V^(r,t)+ffi 



(80) 



Zi = Zi(r,t), 



(71) 



4>{r,t) denoting a suitable scalar potential, so that the 
first two force terms [in Eq. (|79p ] can be represented as 
— Vp + f = — Vp r + La, with 



with r and t spanning respectively the sets (configura- 
tion space) and / (time axis) . In particular the functions 
Zi(r,t) are assumed at least continuous in all points of a 
closed set fi x I (extended configuration space), with 
closure of O. In the remainder we shall require that: 

1. Q is a bounded subset of the Euclidean space E 3 
on M 3 ; 

2. / is identified, when appropriate, either with a 
bounded time interval, i.e., I—]to,ti[ C K, or with 
the real axis R; 

3. in the open set ilxl the functions {Z} are as- 
sumed to be solutions of an appropriate closed set 
of PDE's, denoted as fluid eguations; 

4. by assumption, these equations together with ap- 
propriate initial and boundary conditions are re- 
quired to define a well-posed problem with unique 
strong solution defined everywhere in fix J. 

The deterministic description of an incompressible NS 
fluid is provided by the fluid fields {Z} = {V,p}, with 



p r (r,t)=p(r,t)-(f>(r,t), (81) 

denoting the reduced fluid pressure. Moreover, Eqs. (|72|). 
dZSJ>, LZ4|), ([75j) and (J76]) are respectively the mcompress- 
ibility, isochoricity, NS equations and the initial and 
Dirichlet boundary conditions for {Z} , with {Z a (r)} 
and {Z w (r,t)\ gn } suitably prescribed deterministic ini- 
tial and boundary- value fluid fields, defined respectively 
at the initial time t = t and on the boundary dQ. 



APPENDIX B - TIME EVOLUTION OF THE NS 
FLUID FIELDS 

Let us now show that the phase-space dynamics of 
thermal tracer particles determines uniquely the com- 
plete set of fluid fields {V,p}. To prove the statement, let 
us assume that (at the initial time t a G /) the initial val- 
ues of the fluid fields {V(r o , t Q ),p(r o , t a )} are prescribed 
for all r Q £ f2. The let us require that at any time t > t a 
(with t S J) an arbitrary position r (this includes also the 
case in which r belongs to the boundary <9f2) is reached 
(at least) by two thermal tracer particles with phase- 
space trajectories. Denoting by x^(i) = {r.;(i) = r, Vi(t)} 
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(for i = 1,2) their phase-space states, Eqs.([2])-([5]) then 
imply that for i = 1, 2 it must result 



Vi (t) = V(r, <) + ni (r, t)v th {r, t), 



where the unit vectors rij(r, t) are, by definition, solutions 
of the initial value problem ^ for the initial conditions 

n(r;(i ),t a ) = ni(r oi ,to) and v t h{r,t) = J 2pi p*'^ ■ If we 
require, in particular, that ni(r,t) = — n 2 (r, t) the same 
conclusion is immediate, because then 



V(r,i) 

vth(r,t) 



vi(*) +v 2 (*) 
2 

|vi(t)-v 2 (t)| 



(83) 
(84) 



where the second equation yields p(r, i) in terms of 
Eq. ([24l . Let us now assume, instead, more generally that 
there results ni(r,t) ^ — n 2 (r, t), n 2 (r, t). In this case it 
follows 



V(r,t) 



|viW-V2(*)| 

|ni(r,t) - n 2 (r,f)| 

Vl(t)+V 2 (t) 



[(ni(r,t) + n 2 (r, (r, t). 



Therefore, again, the fluid fields are uniquely deter- 
mined at any (r, t) € £1 X I. 
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